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Abstract 

We discuss the influence of an external field on energy transport in 
superfluid. He-II is not isothermal in presence of Earth gravity; instead, 
it supports finite temperature gradient given by a Fourier-like equation. 
We calculate asymptotic behavior of the effective heat resistance in the 
vicinity of the A-transition. 



1 INTRODUCTION 

Superfluid is known to be (and it was originally referred to as) a super-heat- 
conductor. The heat in superfluid is transfered by the normal flow. This is 
fundamentally different from the Fourier law thermal conduction in ordinary 
liquids. Instead, it has a nature of convection. In uniform superfluid the flow 
(and the heat transfer) is dissipationless, i.e., it happens at constant tempera- 
ture. 

Experimentally observed dissipation in bulk uniform superfluid (i.e., when 
the friction with walls is neglected) is believed to be due to the vortex creation. 
This scenario leads to the temperature gradients with complicated nonlinear 
dependence on the supported heat flux. We show that gravity can also explain 
this dissipation. 

Earth gravity, like any external field, destroys the system homogeneity: it 
produces a pressure gradient across the sample. This inhomogeneity is a cause 
of the normal flow dissipation. In other words, gravity is a friction mechanism 
between the normal flow and Earth. Due to the friction, the liquid holds a finite 
temperature gradient whose leading term is proportional to the energy flux. 
Unlike natural convection in ordinary fluids, this scenario of energy transfer is 
described by a Fourier-like equation. 

To investigate the effect quantitatively one needs to correct superfluid hy- 
drodynamic equations by adding external field. In Sec[21 we write down these 
equations. The form of dissipative terms however is not self-evident. To ascer- 
tain correct appearance of the equations we use Landau-Khalatnikov superfluid 
hydrodynamic theory in a free-fall frame of reference. Translating all terms to 
the laboratory frame we unambiguously get sought equations. 
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In SecO we apply obtained equations to the problem of energy transfer 
in superfluid. For simplicity we restrict ourselves to one-dimensional vertical 
setup where influence of walls is neglected and the energy goes either up or 
down. With the help of linearized hydrodynamic equations we find temperature 
gradient induced by the heat flux. In linear regime the flux and the gradient 
are proportional to each other. We conclude the section by an expression for 
the effective heat resistance of this system. 

Effective heat resistance derived in Sec|3]becomes larger in the vicinity of the 
critical temperature. In SecQJwe show that it diverges at A-point and calculate 
its critical behavior. 

2 HYDRODYNAMICS 

Ordinary superfluid hydrodynamic equations are modified by external fields. To 
find out the effect of gravity we use the equivalence principle of general relativity. 
For the uniform gravity field of Earth we define a non-inertial free-fall frame of 
reference K . It has constant acceleration g with respect to the laboratory frame 
K(t,x l ). We also demand that at t = this frame has zero velocity relative 
to K. This ensures that hydrodynamic variables are equal in two frames. In 
the frame K no gravity is present and (in Newtonian limit) the superfluid is 
governed by conventional equations [2] 
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where tilde attributes quantities to K, p is the fluid density, j is the mass flux, 
v s is the superfluid flow velocity, and E is the energy density. The fluxes in £[J 
are given by 

Il ik =v i J k +v k j i -pv k vi+pS ik +n ik , (2) 
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and 

{ v 2 \ ~ dT 

Q l = U + -f) f + ST < + (/' ~ pv k s ) + (f ~ K)4>+ vfr ik - K—.. (6) 

Here v n is the velocity of the normal flow, p is the pressure, fi is the chemical 
potential, S is the entropy density, T is the temperature, and 77, d, £ 2 i C3j K are 
kinetic coefficients. As usual, we took not all the dissipative terms that are in 
principle possible, but only the largest. In general the dissipation in isotropic 
superfluid is described by 13 kinetic coefficients rather than by 5 (A. Clark 1963, 
see discussion in Ref.pQ, p. 525). 

To rewrite the equations in laboratory frame of reference one must translate 
left-hand side terms in (|TJ to K . We finally get 



dp 
dt 


+ df_ 

dx 1 


= 


Qjk 


dfl lk 


= P9 k 


dt 


dx 1 


dv!: 
~dt ' 


<9$ 

dx k 


= 9 k 


dE 

~bl' 


^d&_ 

dx % 


= J k 9 k 



(7) 



3 EFFECTIVE HEAT RESISTANCE 

For the steady-state flow all time derivatives in J7J) vanish. Consider one- 
dimensional system where all velocities and fluxes are vertical. Let z-axis run 
down along g. In typical heat transfer experiment j^] the energy is carried by 
the superfluid counterflow, where normal and superfluid velocities are directed 
oppositely to each other to keep zero net mass flux. The condition j = is sat- 
isfied along the height of the cell. Under these circumstances the first equation 
in (0 is trivially satisfied and the other three have the form 

n' = pg 

0' = 9 (8) 
Q' = 0, 

where dashes denote the d/dz derivative. 

Within linear (on velocity and temperature gradient) approximation the 
fluxes in JHJ are given by 

U = p + Ci{pv n )' -C2v' n -^r,vL, (9) 
= A* + Ca(p«n) , -Ci<, (10) 
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Q = TSv n — kT . 



(11) 



Keeping lowest nontrivial gradient terms (thus assuming the effect of g to be 
small) and excluding p! by means of thermodynamic identity p dp = dp — S dT, 
we obtain from ©, © and (JTDJ 

A'v' n + Bv u + ST' = 0, 

where we introduced 

A = 2pCi - ~r? - C2 - P 2 C 3 , 

s = (Ci-pC3)p' + (Ci-pC3)p". 

Finally employing JHJ and (| 1 1 1) we obtain a Fourier-like equation 

T' = -qQ, (12) 
where the effective heat resistance g is determined by 

We see that in linear approximation the temperature gradient is propor- 
tional to the heat flux and does not depend on its direction. The effective heat 
resistance is of the second order in system inhomogeneity, i.e., g oc g 2 . If higher 
terms are taken into account g will be different when heated from above or from 
below. 

4 CRITICAL BEHAVIOR 

The effect discussed in previous section will be present in any Earth-based energy 
transport experiment, but it seems to be beyond experimental accuracy unless 
measured in the vicinity of the A-point. To find critical behavior of the effective 
heat resistance we need to estimate kinetic coefficients entering (IB} . 

This can be done using the ideas of dynamic scaling j^. The specific heat 
index for helium is very small 5 , we further let it equal to zero. The superfluid 
density p s follows Josephson scaling relation [H] p s oc r" 1 oc i 2 / 3 , where r c is the 
correlation length and t is the reduced temperature t = T\—T* 

Second sound damping at the boundary of the critical region is inversely 
proportional to the correlation length 

*The same letter t is used to represent the time and the reduced temperature. This should 
not lead to confusion however: the meaning is clear from the context. 



4 



Here p n = p — p s , uj and u 2 are the complex frequency and the velocity of the 
second sound respectively. The latter is determined by the expression 

S dT p s 1/3 

Further assuming that at the critical region boundary all terms in i|14[l make 
contribution of the same order to the second sound damping we get for the 
kinetic coefficients 

77 cx Ci oc C 2 °c C3 °c r c u 2 Ps X oc t~ x . (15) 
To find the main term of the space derivatives in (jlrifl we note that 

Here the derivative dT\/ dp is taken along the phase transition curve. Substi- 
tuting this in and using (|15|l we get^ 

gOC^. (16) 

Having shown that finite temperature gradient in superfluid may exist, it is 
instructive to find the maximum gradient possible. Obtained results on effective 
heat resistance (being derived within linear approximation) correspond to small 
temperature gradient and velocities. They are therefore not applicable to the 
critical velocities region. We will try, however, to use them as an estimate 
beyond the limits. It is shown in Ref. [5j that the critical energy flux Q c in 
the vicinity of the A-transition behaves roug hly as Q c oc i 4 / 3 . Using CH) we 
see that the maximum temperature gradient T c ' increases when the temperature 
approaches T\ as T' c oc < -2 / 3 . 
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